Normal state pseudogap and (tt, 0) feature in the underdoped high-T c cuprates. A 

microscopical theory. 
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We show that a simple 2D electron system on a square lattice with hoping between more than 
nearest neighbors exhibits in the presence of electron spin exchange interaction properties strikingly 
similar to those observed in the underdoped cuprates It is a normal state pseudogap, its behaviour 
with doping and T, a form of the spectrum around (0, tt) etc. The motor is the electronic topological 
transition which due to two- dimensionality gets new features with respect to well-studied 3D case. 
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Many experiments for the hole doped high T c cuprates 
provide an evidence for the existence of a pseudogap in 
the underdoped regime above T c , see for example 
Its value increases with increasing doping distance from 
the optimal doping, 8 opt — 8, the pseudogap exists until 
rather high temperature T* ap (S) which increases in the 
same way with S op t — 8 while the critical temperature 
of superconducting (SC) transition decreases that makes 
doubtful its SC origin. ARPES experiments ||-|5|] show 
that the gap opens for wavevectors around (0,7r), the 
spectrum gets a flat shape and it " stops" at some thresh- 
old wavevector in the direction (0, 7r)-(7r, tt), the spectral 
functions have very damped form for k in the vicinity 
of (0,7r) (so called (0, 7r)-features). The problem of the 
electronic anomalies above T c is considered today as a 
key problem which should help to better understand the 
microscopic origin of high-T c superconductivity. 

Two approaches are possible to attack the problem. 
A phenomenological one which uses an experimental 
information, for example that information (from neu- 
tron scattering and NMR) which tells us that significant 
bosonic field (corresponding to magnetic fluctuations) ex- 
ists above T c just in the underdoped regime. This pushes 
up the idea that the anomalous electron spectrum can be 
a result of fermion-boson interaction. The theories 
are of such type; the bosonic fluctuations are presumed 
to be a memory about AF insulating state existing at low 
doping. There are also theories which consider bosonic 
field related to SC instability H and polaron effect Jl0| ] . 

The second type of approach, microscopical, should 
obtain the bosonic field if it is a fermion-boson interac- 
tion which is at the origin of the phenomena and should 
explain why this bosonic field is strong in the so 
wide range of doping and temperature quite far 
from both AF and SC ordered phases. This way is of 
course preferable as it allows to reveal a nature of effects. 
The theory which we start to develop in jllj is of the 
second type. We have considered a metallic state in 2D 
electron system on a square lattice. We have shown firstly 
that an electronic topological transition (ETT) occuring 
in the system of noninteracting electrons at electron con- 
centration n c when Fermi level (FL) crosses saddle-point 
(SP) is quite untrivial in the case of hoping between more 
than nearest neighbours, t' ^ or/and t" ^ 0,... (more 



general in all cases excluding the perfect nesting). Cor- 
responding quantum critical point (QCP) 8 = 8 C , T = 
(8 = 1 — n) combines two independent aspects of criti- 
cality, the ordinary one related to logarithmic singulari- 
ties in thermodynamical properties and the untrivial one: 
this point is the end point of the critical line related 
to Kohn singularities in 2D electron polarizability [ pd| . 
Quite trivial consequence of this ETT is a developing of 
density wave (DW) and SC instabilities around the point 
8 = 8 C , T — in the presence of interaction and therefore 
an existence of critical bosonic fluctuations in the disor- 
dered metallic state above them. The untrivial aspects 
concerning to DW degrees of freedom are: (i) strong as- 
symetry between regimes 8 < 8 C and 8 > 8 C and (ii) very 
long (in doping and temperature) memory about DW 
instability in the disordered state on one side of ETT, 
8 < 8 C . This encourrages us to continue our analysis and 
to consider now the interaction of bare electrons with 
DW bosonic field that is the aim of the present paper. 

The basic equations are following. The spectral func- 
tion corresponding to the renormalized Green function, 
A(k, ft) = lim 7 _o \ImG((k, ft + i-y)\, is given by 



A(k,Q) = 



|JmS(k,f2)| 



[ft - £k CT - #e£(k, ft)] 2 + [Jm£(k, ft)] 2 
with the irreducible part S(k, iCl n ) determined as 
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[We work with electron spin exchange interaction, the 
same which leads to SC instability [O] and SDW in- 
stabilty [p~T[| l Being practically standard for the case of 
fcrmions interacting with bosonic field, the equations 
lead to results which explicit form crucially depends on 
explicit forms of the bare electron spectrum and of 
the bosonic Green function (or corresponding suscepti- 
bility %). In our approach we use the susceptibility cal- 
culated for the considered 2D electron system in the dis- 
ordered state above both SDW and SC instabilities. And 
we use the spectrum which posseses SP : in the form 
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£ko- = —2t(cosk x + cosky) — At cosk x cosk y (3) 
for numerical calculations and in the form 



£k<r = fiktr - A* = _ ^ + _ bk 



(4) 



fc y are distances 



for analytical estimations. In (^), k 
from SP wavevector (0, tt), a = t- 2t' , b = t + 2t' . The 
important parameter is defined as Z — /i — e s (e s = 
—At' ft is SP energy) and measured the energy distance 
from ETT : Z oc 5 C — 5 |n| . This parameter determines 
a new energy scale in the system. As we will see, 
results for £ and therefore for the renormalized electron 
Green function are very different from those obtained in 
Jg,[2/ where the SP in the bare spectrum was ignored. 

As shown in pH , the susceptibility in the regime 
6 < Sdw{T = 0) of the disordered metallic state can 
be roughly presented in the form 



x(q,w) 



AJ ( k 2 



(1 - K 2 ) + iC(uj)u 



A(q- 
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valid for low energies to <C lo c oc 5 C — 5, with C(0) oc l/Z , 
the parameter k 2 describes a proximity to the DW in- 
stability line. The important parameter is the energy 
for which Jm^(q,u) is maximum, the param- 
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C(0) 



eter which determines the second energy scale in the 
system. The behaviour of n 2 , C(0), and ujo, as functions 
of S and T is quite untrivial, we discuss this point later 
on. We use the form (^|) for analytical estimations while 
RPA form with the electron-hole loop calculated based 
on the full spectrum (||) jll| for numerical calculations. 

Analytical calculations for T=0 performed with the 
spectrum (^) and susceptibility for k=ksp=(0, tt) 
give for SI corresponding to the bare spectrum ek SJ3 = — Z 



ImT,(k S p, SI = 6k S p) 00 



C(0) 



(to = (1/a + l/b)/2) while for low energies SI <C wo they 
give : /TOE(k, SI) oc SI 2 . Numerical calculations per- 
formed with the full spectrum (||) and with the calculated 
RPA susceptibility give the results shown in Fig.^Ja. Both 
numerical and analytical calculations show that imE 
taken for fixed Z becomes singular at SI = ek SP = — Z 
when k 2 — > and that for very low energies the Fermi 
liquid behaviour survives. The corresponding behaviour 
of ReT, is presented in Fig.|l|b together with the line SI — 
since the equation SI — e^ — i?eS(k, SI) = determines new 
poles. One can see that there are three poles for small 
k 2 . The position of the pole at the intermediate energy 
is close to the position of the pole in the bare spectrum 
while the two external poles are precursors of the new 
bands in the ordered DW phase, see jl3). The impor- 
tant point is that the spectral function exhibits maxima 
approximately at the energies corresponding to the new 
poles while it is minimum at the position of the 
bare pole, see Fig.|c. This occurs due to the singularity 
(or maximum in the case of finite k 2 ) in |JmS| just at SI 



corresponding to the position of the bare spectrum. Such 
a behaviour means an appearence of two new modes of 
the spectrum instead of one in the bare spectrum. For in- 
termediate k 2 the poles at negative SI disappear, however 
the peak in the spectral function at negative SI survives 
being now a resonance peak. 
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FIG. 1. /mS, i?eS and spectral functions as functions of 
D. for k = k S p and different values of k 2 = 0.05,0.01,0.001 
depending on t/J. (T = 0, Z/t = 0.2, t' ft = -0.3.) 

Calculations of 7ro£ for k along different directions 
4> traversing SP show that the effect of the maximum 
in |/toS| is preserved for wavevectors located in some 
distance Ak((j>) from (tt, 0), see for example Fig.|| corre- 
sponding to the direction tt(1 — k, k) which traverses the 
"hot spot" . The SI 2 behaviour at low SI is also preserved. 
[One should note that there is no square-root behaviour 
|| for SI > lu in the presence of SP.] This picture is 
valid for all directions including (0, 7r) - (tt,tt) and (0, 7r) 
- (0,0). The maximum of |7mE| (when exists) occurs at 
SI* = —\Z + auiQ + b(Ak((b)) 2 }. Above some threshold 
value of Ak*(<f)), ImT, becomes almost unstructured. 
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FIG. 2. /mS(k, Q.) as functions of £1 for wavevectors in 
the direction tt(1 - k, k). T=0, Z/t=0.3, t/J=1.8 (k 2 =0.066, 
uo/t= 0.08 

There are number of important consequences of such 
a behaviour. The first is the form of the spectrum in 
the vicinity of (tt, 0) : the bare spectrum splits into two 
branches, and a gap opens, see Fig.[|a for the spectrum 
obtained from maxima of spectral functions. Secondly, 
the spliting into two branches and the gap disap- 
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pear above some threshold wavevector when one goes 
far away from SP, the bare spectrum is restored. The 
corresponding density of states (DOS) calculated for the 
same conditions as the spectrum in FigJ^a is shown in 
Fig.^b. It exhibits a pseudogap and two new peaks in- 
stead of the SP peak in the bare spectrum. 
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FIG. 3. Electron spectrum along directions (0, 7r) — (tt, tt) 
and (0, tt) — (0, 0) in two temperature regimes, I (T<&jo, 
u>o<^l) and II (T > uq) for fixed Z. Calculations are done 
for (a) T=0 and (c) T/t=0.16, in both cases for Z/t=0.3 
(8 = 0.1). For the regime I the corresponding DOS is shown 
for comparison (b). Results for the regime II are compared 
with ARPES data [3] for the underdoped BSCO (d). The 
dashed lines in (a), (b), (c) correspond to the bare spectrum. 

Thirdly, being directly related to the energy scale 
Zoc5 c — S, the gap increases with Z or by other words 
with decreasing doping, see Fig^J. The coefficient of 
proportionality depends on interaction increasing with 
J/t. Fourthly, due to high value of Im£| for corre- 
sponding to the lower branch this branch is strongly 
damped in the vicinity of (0, tt). The damping decreases 
when one goes away from SP and the well-defined bare 
spectrum is restored. All these features are in a very 
good agreement with ARPES. There is, however one im- 
portant difference: the upper well-defined branch charac- 
terized by new SP located almost on FL is not observed 
experimentally. This could come from finite T effect. 




FIG. 4. T = pseudogap (determined at (0, tt)) as a 
function of doping for different values of t/J ( t' /t = —0.3). 
Each curve starts at 8 = Sow which value is determined by 
t/J. 

To check let's study a behaviour at finite T. After the 
analysis similar to that performed for T = 0, one can real- 
ize that there are four different regimes in T. For the first 
three, the lower branch in the spectrum keeps its main 
features while the upper branch changes its behaviour. 
Namely, in the low temperature regime I, T « w , 
k 2 « 1, the energy cb is negative and the upper branch 
traverses FL at some point in the direction (0, tt) — (tt, it) 
as in Fig.||a. FS exists being open but quite close to 
the critical form. In the crossover regime, T ~ ujq, one 
has cb ~ 0, the new SP in the upper branch traverses 
FL while this branch becomes much more damped than 
in the regime I. In the regime II, T > uq, the SP energy 
moves to positive fi, eg > 0, and what is most important 
the upper branch does not traverses FL anymore 
(see spectral functions in Fig.^Ja, 0b and the spectrum 
obtained from their maxima in Fig!yc) : Fermi surface 
disappears. And finally in the regime III, T < ujq while 
both parameters are not very small, the spectral func- 
tions have a form shown in Fig.^c,d. The gap is filled up, 
the bare spectrum is restored although it remains quite 
damped at (tt, 0). 
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FIG. 5. Spectral functions calculated for two directions 
traversing SP (n, 0)-(tt, tt) and (-7T, 0)-(0, 0) and for two tem- 
perature regimes, II and III for fixed doping 5=0.1 (Z/t=0.3). 
The wavevectors are done in tt units, t' /t=— 0.3, t/J=1.8 

Areas in T ~ S plane where each regime takes place de- 
pend on behaviour of uq as a function of S and T. This 
behaviour is analyzed in fO]. We remind that w (taken 
at fixed doping S < &dw)i slightly decreases with increas- 
ing T below Tft e (5) = a(5 c — 5), remains constant until 
another temperature, Tj m (S) — (3(S C — S) and then in- 
creases with T. This leads to the phase diagram in Figjfj. 
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where we show one more regime, the regime IV, T ujq. 
For it a gap opens for all k in the direction (1 — k, k), 
the spectrum has the form as in the ordered phase, see 
Fig. 6 in ju|. The origin is a peak in |/mS(k, 0)| at f2= 
£k+Q AF existing for all k. This regime corresponds to the 
quasistatic regime analysed in (although the authors 
did not note that the gap opens for all wavevectors along 
(k, 1 — k)). It is important to emphasize that it takes 
place only in the intime vicinity of AF instability 
whatever its origin (SDW or the spin-localized as in M). 
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FIG. 6. Calculated phase diagram for t — t electron sys- 
tem in the presence of J-interaction. Tow is the line of DW 
instability [11], T sc of SC instability [12], different crossover 
regimes are discussed in the text, t /t = —0.3, t/J = 1.8. 

All regimes discussed have been obtained in ignoring 
a superconductivity. In the presence of the latter (SC 
instability also develops around the QCP as discussed 
in |lj,|ll[]) the regimes I, IV are unacessible being cov- 
ered by SC phase. It is the regime II which occupies the 
most interesting part of T — 5 plane and which should 
be compared with experiment for the pseudogap regime. 
The corresponding ARPES spectrum for the underdoped 
BSCO H is shown in Fig.||d. The agreement is excellent. 



untrivial f2 dependence, for two wavevectors located 
not far from (n, 0) in two directions (tt, 0) — (n, n) and 
(tt, 0) — (0,0) to emphasize the absence of crossing of 
FL, and for the wavevector close to kp in (1,1) direc- 
tion to emphasize the ordinary FL form of the response 
function. The behaviour is very close to that observed 
experimentally ||. All these regimes are related to the 
DW fluctuations. The question arises about a precursor 
of the superconductivity. Do corresponding fluctuations 
succeed in opening a pseudogap. The answer is that it is 
indeed so but only in the intime vicinity of T SC {S) in the 
same way as for the regime IV with respect to the line 
Tdw{5); the details will be discussed elsewhere. 

Summarizing, we have found that a simple 2D electron 
system on a square lattice exhibits in the presence of elec- 
tron spin exchange interaction properties strikingly sim- 
ilar to those observed in the underdoped cuprates, the 
normal state pseudogap, its behaviour with doping, its 
existence below T* increasing with decreasing <5, the 
shape of the spectrum around (0, ir) etc. : the fact which 
is quite impressive in a view that we did not use any ex- 
ternal hypothesis or adjastable parameters. The motor is 
ETT which due to two-dimensionality gets new features 
pd| with respect to well-studied 3D case, sec [ pT|Jl5| ] . The 
effects disappear in the nesting case, t' = t" = ... = 0: the 
regimes 6 < 5 C and 6 > 6 C become symmetrical, for both 
regimes the bosonic field is incommensurate and rapidly 
weakening with increasing \S C — S\ and T jyj. The re- 
sults are in a good agreement with exact diagonalization 
study for t — t' — J model (5). The general picture is 
not sensitive to values of microscopical parameters, t' /t 
and t/J (only positions of the characteristic points, 5 C , 
6p W , 5~jj W depend on them). The theory predicts a quite 
interesting behaviour for the spectrum and spectral func- 
tions for positive uj that would be interesting to check by 
photoemission experiment. 
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FIG. 7. The function A(k, lu)u f (lu) for Z/i=0.3, T/t=0.16 
and four characteristic wavevectors. t'/t=— 0.3, t/J=1.8 

With the same aim to compare with experiment we 
show in Fig.0 the response function A(h, oj)n F (uj) cor- 
responding to that measured by ARPES for four most 
characteristic wavevectors : for (tt, 0) to emphasize the 
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